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Abstract 



We introduce a general approximation scheme in order to calculate gauge invariant observables in 
the canonical formulation of general relativity. Using this scheme we will show how the observables 
and the dynamics of field theories on a fixed background or equivalently the observables of the 
linearized theory can be understood as an approximation to the observables in full general relativity. 



Gauge invariant corrections can be calculated up to an arbitrary high order and we will explicitly 
^ , calculate the first non-trivial correction. 

Furthermore we will make a first investigation into the Poisson algebra between observables corre- 
' sponding to fields at different space-time points and consider the locality properties of the observables. 

o' 

q-i 1 Introduction 



One of the most important issues for a quantum theory of gravity is the construction and interpretation of 
the observables of the theory. We will address this issue in the framework of the canonical theory, where 
the gauge independent observables are called Dirac observables. So far - apart from the ADM charges 
PP in the context of asymptotically flat space-times - no Dirac observables are known explicitly for pure 
gravity. Moreover in [2] it is shown that with the exception of the ADM charges, Dirac observables have 
to include an infinite number of spatial derivatives. Hence we expect that it will be tremendously difficult 
to calculate Dirac observables and that the only resort may be to develop approximation methods for 
Dirac observables in order to make physical predictions. One proposal [3] for an approximation scheme 
is an expansion in the inverse cosmological constant. 

In this work we will propose a general approximation scheme for Dirac observables, based on the 
concept of complete observables JI . A complete observable ^[/T]( r ) is a special kind of relational 
observable: Using some of the dynamical variables T as clocks, the complete observable gives the value 
of some other dynamical variable / at that instant at which the clock variables assume certain values r. 

The complete observables can be computed using a power series [SUS] in the (phase space dependent) 
clock variables. This power series is a natural starting point for a perturbative approach and indeed we 
will use it to obtain an approximation scheme by expanding the series in the "fluctuations" around some 
fixed phase space point. 

We will apply this approximation scheme to general relativity. In doing so, we have to make certain 
choices - the most important being the choice of the clock variables. Here our guide line is that we want 
to have a good approximation to the observables of field theory on a fixed background which in this work 
will be the flat Minkowski background. As we will see this results in observables which in the zero gravity 
limit (i.e. for k = and vanishing gravitational fluctuations) coincide with the usual observables of field 
theory on a flat background. The gravity corrections can be calculated explicitly order by order and are 
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connected to the standard perturbation theory. The first order observables are given by the observables of 
linearized general relativity, hence this method gives us a precise understanding of the observables of the 
linearized theory, for instance the graviton, as approximations to observables of the full theory. Moreover 
the approximation scheme in this work gives a precise proposal how to compute higher order corrections 
to the observables of the linearized theory. These higher order corrections are in a consistent way gauge 
invariant to a certain order - to make these corrections completely gauge invariant one would have to 
add terms which are of higher order than the corrections themselves. One important point is that we do 
not only manage to approximate well the fields at one specific time, but obtain also an approximation to 
the fields at arbitrary times, which makes the extraction of dynamical information much easier. 

This brings us to the second issue mentioned above, namely the interpretation of the observables, 
in particular in view of the quantum theory. One question for instance is, whether it is possible at all 
to construct Dirac observables which give the standard (local) field observables in some limit and what 
kind of locality properties these observables have. Having a (perturbative) computation scheme for Dirac 
observables at hand allows one to examine these issues in more detail. 

Deviations from the standard field observables could result in fundamental uncertainties for the ob- 
servables in a quantum theory including gravity, as is argued in [J] using relational observables. One 
crucial aspect for these deviations is that if one uses relational observables one has to specify some of 
the dynamical fields as clock variables. However these clock variables are dynamical. Since measurement 
involves always a disturbance of the system this means that one has also to expect a disturbance of the 
clock variables if one would measure a complete observable. Another aspect is, that if we use matter 
fields as clocks these matter fields are coupled to the gravitational field and therefore influence all the 
other fields. In the classical theory these influences will show up for example in the Poisson algebra of 
Dirac observables corresponding to measurements at different space-time points. Here it is important 
to consider space-time points at different times, since equal-time Poisson brackets will be zero for the 
cases we are interested in. The Poisson algebra will be reflected in uncertainty relations in the quantum 
theory for these measurements. Connected to this is the question whether this Poisson algebra satisfies 
any locality properties, for instance whether fields at spatially separated points Poisson commute. In this 
work we will make a first investigation into the classical Poisson algebra of observables corresponding to 
measurements at different space-time points. 

However the Poisson algebra and the uncertainty relations which follow from this algebra will depend 
on the choice of clock variables and it is important to understand in which ways this choice matters. Let 
us consider a very simple example, namely two parametrized particles with the Hamiltonian constraint 

C= Pt + ^- + ^- (1.1) 

where pt is the conjugated momentum to the time variable t and pi,p 2 are conjugated to the two position 
variables (71,(72- A natural choice for a clock variable is t and we can ask for the position of the first 
particle at that moment at which t assumes the value r. We will denote this observable by Fi qi . t -\(r) and 
it can be easily calculated to be 

F [qi;t] (T) = qi + ^(T-t) . (1.2) 

TO 

It Poisson commutes with the constraint and is therefore a Dirac observable. The Poisson bracket of two 
observables (n) and i 7, [ gi;t ](r 2 ) at two different clock values t\ and r 2 is phase space independent 

{F [qi , t] (T 1 ),F [qi , t] (T 2 )} = ^-(r 2 -T 1 ) . (1.3) 

Now one could also choose the position of the second particle as a clock variable and ask for the position 
of the first particle at that moment at which the second particle has position r'. The corresponding 
observable is 

F [qi ; q2] (r') = qi + ^ — (r'-q 2 ) . (1.4) 
mi P2 

If one ignores that r and r' refer to different clocks, (|1.4fl looks of course quite different from (|1.3|) . 
However if one takes into account that the value r is reached at that moment at which 

r' = F [q2 . t] =q 2 + ^-(r-t) (1.5) 
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and uses this to replace r' in (|1.4|) one will get back to l|1.3|) . In so far both choices of clock variables give 
us the same time evolution if one takes into account the "translation" (|1.5|l between the clock readings t 
and t'. This is only a classical consideration, the quantization of i|I.4[l will have quite different properties 
from the standard position operator due to the difficulties involved in quantizing the inverse momentum 
p-2 in (|1.4(l . see [El El and references therein. 

However the Poisson bracket between two observables with respect to the clock variable q^ is given 

by 

{F\^M^ F \^M)}= [— (T2-r[)][^][l + ^-^] ■ (1.6) 
11 L mi P2 mi P2 

The first factor in square brackets on the right hand side of ljl.6|l is equal to our previous result (|1 , the 
second factor is due to the translation between r- and r'-parameters. The third factor can be seen as 
a correction term resulting from the use of an unusual clock variable. The correction is proportional to 
the kinetic energy of the particle observed divided by the kinetic energy of the particle used as a clock. 
In the limit of large energy of the particle used as clock compared to the observed particle we get back 
to the previous result (|I.3(I . This corresponds to the intuition that the second particle is a good clock if 
its velocity is large, i.e. if it has large energy. Interestingly, if one attempts to quantize operators of the 
kind (|1.4I) one derives an inherent uncertainty inverse to the kinetic energy of the particle used as clock 
U3- To take the limit of large energy for the clock variables is problematic in general relativity because 
of backreaction and ultimately black hole formation. 

We see that different choices for the clock variables lead to different predictions for the uncertainties 
for the complete observables. For general relativity the question arises what kind of clock variables are 
available. For matter clocks and clocks built from curvature scalars we expect a behaviour similar to 
()I.6(I . These clocks have to have a certain energy in order to be good clocks and therefore it is in question 
whether the limit to flat space can be performed. The question arises whether there are any "good" 
clocks near flat space built from the gravitational degrees of freedom and what properties the Poisson 
algebra of the corresponding complete observables has. 

1.1 Overview 

In section|21we will explain how given a first class constraint set, we can construct new constraints with a 
different constraint algebra: The Poisson brackets between these new constraints are at least quadratic in 
the constraints. These new constraints are needed for an explicit power series for the complete observables. 
We will also introduce an iterative method to obtain a constraint set which is Abelian. These Abelian 
constraints will simplify calculations for the case of general relativity and facilitate the interpretation of 
the resulting formula. 

The power series for complete observables serves as a basis for our approximation scheme, which we 
will define in section [3J The basic idea is to keep in the power series only terms up to a certain order 
q in the fluctuations around a certain "background" phase space point. For special values of the clock 
parameters r this will result in only finitely many terms in the power series. 

In the next section 0] we summarize some material necessary to apply this approximation scheme to 
general relativity and in particular we define our clock variables. In[3]we will find that this approximation 
scheme also allows one to compute complete observables corresponding to a one-parameter family of 
clock parameters, that is to different times. This facilitates the dynamical interpretation of the complete 
observables. 

Next we will consider the second order approximation in El and give explicit formulae for the terms 
one has to calculate. These calculations are performed for a Klein-Gordon field coupled to gravity in[7| 
The Poisson algebra between the resulting approximate observables is considered in section |S] 

2 Perturbative Abelianization of the constraints 

In this section we will collect some facts about first class constrained systems that will be important later 
on. More concretely we will show how one can redefine the constraints such that these new constraints 
have particular nice properties, which are key for the construction of the complete observables. In the 
second part we will introduce a method to obtain an Abelian set of constraints. 
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To simplify notation we will restrict ourselves to a finite dimensional phase space, however the results 
are easily generalizable to an infinite dimensional phase space. 

Let {Cj}JL 1 be a set of independent 1 first class constraints on a phase space described by canonical 
coordinates {X a }^ l , l 1 . Choose a phase space point Xq := {Xq}^L 1 satisfying the constraints and intro- 
duce new canonical variables as "fluctuations" x a := X a — Xq around the phase space point. Rewrite 
the constraints in this new variables and define the first order constraints ^'Cj as the part of the con- 
straints linear in the fluctuation variables. We have to assume that these linearized constraints are also 
independent from each other, in particular they should not vanish. 

Now the linearized constraints are Abelian ^U] 

{«C,-,«C fc } = (2.1) 
as can be seen by expanding the full constraint algebra 

{ V) Cj + ^C 3 +..., Wft + ^C k + ...}= (Wf jk + + . . .){WQ + + ...) (2.2) 

order by order. Here the superscript ^ left to a symbol denotes its m-th order. Hence one can find a 
set of clock variables {T K } r ^ =1 canonically conjugated to the linearized constraints, i.e. such that 

{T K ,WCj} = 6? . (2.3) 

This set of variables can be extended to a canonical coordinate system 2 by finding an additional set of first 
order quantities {Q }2=j™ an d {Pd^ZT having vanishing Poisson bracket with the linearized constraints 
and with the clock variables and satisfying {Q d , P e } = 
Define the phase space dependent matrix 

Af := {T K , Cj} = Sf + Bf (2.4) 

where B^ is at least a first order quantity. Because the zeroth order of is given by the identity matrix 
we can invert Aj order by order. We will assume that A^ is invertible in some phase phase space region 
around Xq, i.e. that the series defining the perturbative inverse of A^ converges in this region. In the 
following we will restrict our considerations to this region of the phase space. 

Using the inverse of Aj we define a new set of constraints (equivalent to the old one at all phase 
space points where the inverse can be defined) by 

C K = (A-%C J . (2.5) 

This new set has the following properties: First it is easy to see that 

{T K , Cm} = S M + \ M N C N ~ S M (2.6) 

where Af^ = {T K , {A^f^Af are a set of phase space functions and the sign ~ means weakly equal, 
that is equal up to terms vanishing on the constraint hypersurface. 

From equation l|2.6|l one can prove that the new constraints are "weakly abelian" , i.e. they Poisson 
commute up to terms quadratic in the constraints: Compute {{T K ,Cm},Cn} first directly and then 
using the Jacobi identity. Comparing the two results one can conclude that the structure functions fjfj 
defined by {Cjc , Cj} = f¥j^M have to vanish on the constraint hypersurface. 

Property (|2.6(l of the constraints Ck will be the key property in order to be able to write down a power 
series for the complete observables. In the following we will construct a set of constraints Ck, which 
satisfies {T K , Cm} — 5 M exactly, i.e. also away from the constraint hypersurface. These constraints will 
simplify very much the calculations in the later sections. 

Note that the flow generated on the constraint hypersurface by a constraint Cl does not change if we 
add to Cl a combination at least quadratic in the constraints 

{/, C L } ~ {/, C L + fi KM C K C M } (2.7) 

where \i K M is an arbitrary smooth phase space function. 

M.e. the equations Cj = should be independent from each other 

2 In general this canonical coordinate system is only a local one, i.e. it can be only denned in some region around Xq. 
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In particular if we add to the constraints Cx a combination quadratic or of higher order in the 
constraints, equation (|2.(i|) will still hold, the only difference is the exact expression for the phase space 
functions X M . However the exact expressions for these coefficients are not important for the proof that 
the constraints Cl are weakly abelian. 

This can be used in an iterative process to find a set of Abelian constraints (if this process converges) . 
This process can be applied to an arbitrary set of clock variables, as long as the clock variables are 
Abelian. In the following we will explain the iteration step. 

Assume that one has a set of independent first class constraints {Ck }k=i ana ^ a se ^ °f independent 
phase space functions {T K }^ =1 with the property 

{T K ,C M } = 6 M + X M L ^C Ll ---C u . (2.8) 

for some r > 1 and such that 

{C M , C n } = ftf N Lr+1 C Ll ■ ■ ■ C Lr+1 = 0(C('' +1 )) (2.9) 

for some phase space functions fjJ~N r+1 ■ Here (^(C 1 ) means a smooth function g vanishing on the 
constraint hypersurface. Such a function g can be written as a combination of the constraints g = v K Ck 
with some smooth phase space functions Vk, see JTTJ. With 0(C r ) we denote a phase space function g 
which can be written as g = v Ki ' Kt Ck 1 • ■ ■ Ck t with some smooth phase space functions v Kl ' Kr . 

The X M Ll "' Lr in (|2.8|l can be a set of arbitrary (smooth in the phase space region of interest) phase 
space functions, without loss of generality we can assume that they are invariant under permutation of the 
{Li, - • ■ , L,.} -indices. Furthermore we will assume that the \T k Yk=\ are Abelian, i.e. have vanishing 
Poisson brackets among themselves. 

With these assumptions we will show that one can define a new set of constraints 

C K := C K + &'" Lr+l C Ll ■ ■ ■ C Lr+1 (2.10) 
with fx^ 1 Lr+1 symmetric in the {Li, • • • , L r }- indices such that 

{T K ,C M }=S M +X^- LrLr+1 C Ll ---C Lr+1 . (2.11) 
Furthermore these new constraints satisfy 

{C K , Cm} = f KM ' Lr+2 C Ll ■ ■ ■ C Lt+2 = O(C^) (2.12) 

for some phase space functions ] K \i Lr+2 ■ 

Proof: First note that we can find the coefficients X M Ll Lr in 12. 8[) at least up to terms proportional 
to the constraints by taking iterated Poisson brackets between the clock variables and both sides of 
equation l|2.8|l : 



\ M Ll '- Lr = ^{T L \{T L ',{--- ,{T L ",{T K ,C M }---} + 0(C) . (2.13) 



Hence the X M Ll ' Lr are invariant up to 0(C)-terms under permutation of the {K, Li, • • • , L r }-indices. 
(We have assumed that the clock variables are Abelian.) 

Using equation l|2.10|l as an ansatz in equation l|2.11|l and exploiting the assumption (|2.8|l we arrive 

at 

{T K , C M } = S M + X M L ^C Ll ■ ■ ■ C Lr + (r + l)fi M L ^C Ll ■ ■ ■ C Lr + 0(C^) . (2.14) 
Therefore we have to choose 

% 1 '" £r+1 * - A M'" i(r+1 ' - ~^TT)^ TLl ^ TL2 ^- ■ ■ > {T L '\{T L "+\C M } ■■■} (2.15) 

in order to satisfy equation Q2.11[) . 
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It remains to show that the new constraints commute up to 0(C^ r+2 ^ )-terms. To this end consider 
the Jacobi identity 

{{T k ,Cm},Cn}-{{T k ,C n } ) C m } = {{C n ,C m },T k } . (2.16) 

Because of property |2.11l) the left hand side is of order 0(C^ r+1 ^). Using the definition p.lOjl l for the 
Ck in terms of the Ck and assumption (|2.9(l for the Poisson bracket between the constraints Ck the 
right hand side can be written as 

{{C N , C m }, T k } = {gfa Lr+1 6 Ll ■ ■ ■ C Lr+1 , T K } + 0(C^) 

= (r + l)g N ^ K C Ll ---C Lr +0(C^) (2.17) 

where g N \j Lr+1 is some set of phase space functions symmetric in the {L\ ■ ■ ■ L r+ i}-indices. 

Reinserting this result into equation H2.16|l we see that all terms are of the order 0(C^ r+1 - ) ) except 
for (r + 1) g N 1 M Lrh Cl 1 ■ ■ ■ Cl t - Applying r times the Poisson bracket with the clock variables T Kj with 
both sides of the equation, we conclude that the g N * M r+1 have to vanish weakly. This shows that the 
Poisson bracket between the constraints Ck is of order 0(C^ r+2 - ) ). Hence we are able to define a set of 
constraints that have vanishing Poisson brackets up to terms of arbitrary high order in the constraints. 
If the iteration procedure converges it will result in a set of Abelian constraints. 



This method is applicable to an arbitrary choice of Abelian clock variables. Here we are interested in 
the case where the clock variables are canonically conjugated to the linearized constraints. Let {CkYk=i 
be a constraint set such that 

{T K ,C M } = 5^ + 0(C r ) . (2.18) 

Note that during the procedure described above to find such constraints the first order terms of the 
constraints do not change, i.e. ^'CW = ^Cjd 3 M . We have that 

{T K ,C M - m C M }=0(C r ) (2.19) 

is at least of order r in the fluctuations. Therefore if we would rewrite the constraints Cm using the 
canonical variables (T K , ^'Cj, Q d , P e ) introduced above, all the terms of order higher than one and lower 
than (r + 1) are independent of the momenta {^Cj}™ =1 . Also it is easy to see, that the terms of order 
lower than (r + 1) do not change if one iterates the above procedure to get a set of constraints {C'm}m=i 
satisfying {T K ,C' M } = S M + 0(C S ) with s > r. This again allows us to calculate the completely 
Abclianized constraints order by order. 

If the procedure converges for r — > oo we will end up with a deparametrized form of the constraints, 
i.e. we have constraints °°Ck with 

^Ck = {1) C,5 K + E K (T M , Q d , P e ) (2.20) 

where the phase space functions Ek are independent from the momenta ^Cj. One can show |11| . that 
constraints of such a form are Abelian. The Abclianization procedure introduced here allows one to 
compute the deparametrized form of the constraints order by order and is also applicable if the clock 
variables are quite complicated functions. 

In the following we will assume that the Abelianization procedure indeed converges in a finite phase 
space region around Xo and we will denote by Ck the Abelianized constraints. Because of the consider- 
ations above we need only finitely many steps to calculate a finite order of Ck- 



3 Approximate Dirac Observables 

Property (|2.6() of the constraints Ck enables us to construct Dirac observables F[j ]T ](t) for a set of 
parameters {t k }k =1 via a (formal) power series: 

%;T](r) = h { ' ' ' {/ ' ' ' ' ' } ' {rKl ' TKl) ■ ■ ■ {tK " - ^ ■ (3 ' 1} 

r=0 
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It is straightforward to check that the Poisson bracket of this Dirac observable with the constraints 
Ck vanishes at least weakly. The constraints Ck and Ck generate the same flow on the constraint 
hypersurface, therefore one can replace any of the constraints Cm in with the constraints Cm- 

The observable F\f-T]{ T ) restricted to the hypersurface {T K = t k Yk=\ i s equal to the function /. 
This provides the interpretation for F\r.m{j)\ Since it is gauge invariant, one can see it as a gauge 
invariant extension of /. Hence F\f.-n{r) and Fui-tt,(t) coincide at least weakly if / and /' coincide on 
the hypersurface {T K = t k }^ =1 . 

If the constraints generate also time evolution (i.e. if time evolution is a gauge transformation) 
F[f;T]( T i x ) is known as the complete observable associated to the partial observable (and the clock 
variables T K ). Frj. T i(r, x) gives the value of the partial observable / at the "moment" at which the clock 
variables T K assume the values r . 

Neglecting in the power series (|3.1f) all terms of order higher than q in the differences (t k — T K ), i.e. 
truncating the series to its first (q + 1) summands, we obtain a phase space function that has vanishing 
Poisson brackets with the constraints modulo terms of order q in (r — T ) and modulo constraints. In 
this sense we can obtain approximate Dirac observables, the approximation being good near the gauge 
fixing hypersurface {T K = t k Yk=\- However in many situations it may be quite involved to find the 
inverse of the matrix Aj , which is necessary in order to obtain the new constraints Ck or Ck- 

We can also choose to obtain an approximation around a certain phase space point Xo by expanding 
the complete observable in the fluctuation variables introduced above and by using the specific t k 
parameters given by t k = T K (Xq). As already explained this allows us to find the new constraints Ck 
or Ck order by order since we can invert the matrix order by order. 3 Denote by ' 9 '-F]/ ; t] ( t = T(Xq),x) 
the complete observable i*[/ ; T] (r = T(Xq), x) with terms of order higher than q in the fluctuation variables 
neglected. For the calculation of ' 9 '-F[/ ; t] ( r = T(Xq),x) we will need at most the first (q + 1) terms in 
the series If) and the constraints Ck to the q-th order. 

The truncated complete observable ^Fy. T ^(r = T(X$)) will commute with the constraints modulo 
terms of order q in the fluctuations (and modulo terms vanishing on the constraint hypersurface). Hence 
we can call ^F^. t ^(t = T(Xq),x) an approximate Dirac observable. 

If the power series f|3 . 1 1> for the complete observable converges it defines an exact Dirac observable 
which coincides with the approximate Dirac observable ^F^-x] (t = T(Xq)) modulo terms of order (q+1). 
If the power series does not converge in some phase space region, this will be due to the fact that the 
clock variables do not provide a good parametrization of the gauge orbits in this phase space region 
In this case one can try to find a set of new clock variables T' K with a better behaviour in this respect. 
Assume that the complete observable F\f-T'] (t' = T'(Xq)) associated to these new clock variables and 
the partial observable /' := ^F [f . T] (r = T(X )) can be defined. This complete observable will also 
coincide with ^i^.yj (r = T(Xq)) modulo terms of order (q + 1), as can be seen by examining the power 
series (|3.1JI for a complete observable and using that /' Poisson commutes modulo terms of order q with 
the constraints. 



4 Application to General Relativity 

In this section we will collect all definitions in order to be able to calculate approximate Dirac observables 
for General Relativity. We will work with the (complex) Ashtekar variables [2] (and use the conventions 
in |13p. because the constraints are polynomial in these variables. However the formalism is independent 
from the choice of variables. 

The canonical variables are fields on a spatial manifold £ the coordinates of which we will denote 
by {cr a }a =1 . We will assume that E is diffeomorphic to K 3 . The configuration variables are given by 
a complex connection {A J a }^ a=1 where latin letters from the beginning of the alphabet denote spatial 
indices and from the middle of the alphabet sw(2)-algebra indices: 

Ai=T{+(3Ki (4.1) 

where fj — i/2 is the Immirzi parameter, T J a is the spin connection for the triads e J a and K 3 a is the extrinsic 
curvature. The spatial metric can be calculated from the triads by q a b — s J a e^5jk- The conjugated 

3 This applies also to a more general choice of clock variables than the above one, where the clock variables are conjugated 
to the linearized constraints. A sufficient condition is that the zeroth order of the matrix A 1 ^ should be invertible. 
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momenta are constructed out of the triads 

E? = p-h aa ^e jjlh ei\eil (4.2) 

where e aai ° 2 and Cjj 1 j 2 are totally anti-symmetric tensors with £123 = e 123 = 1. The Poisson brackets 
between these phase space variables are given by 

{Ai(a),E b k (a')} = KSiS b J(a,a') (4.3) 

where k = 8ttGn/c 3 is the gravitational coupling constant. 

The constraints are given by the Gauss constraints Gj (a) , the vector constraints V a (a) and the scalar 
constraints C(cr): 

Gj = K-\d a E? + e M A k a E?) 

V a = K- 1 Fi h E] = K-\d a Ai-d b Ai + z jhl A k a A\)E h j 

C = K-^F^ejuEiE^ . (4.4) 

Note that we use the scalar constraint with density weight two here. 

In this work we will choose as our background phase space point Minkowski space Xq :— (A a 3 ,~E b k) 
with 

Aj = , V\ = p- 1 6l . (4.5) 

We will denote the fluctuations around this background by lower case letters: 

aJ=Aj-Aj , e b k =E b k -E b k . (4.6) 

Contracting these fluctuation variables with the background triad or its inverse, we can convert the 
internal Lie algebra indices to spatial indices: 

a a 6 := a/E 6 fe , e c d = e c k E k d . (4.7) 

Since the background metric is flat we can freely raise or lower the spatial indices. The constraints in 
these variables are given by 

G b :=G^ = K ~\d a e a b +/3e bde a ed + /3e bdf e af a a d ) (4.8) 

V a = K-\d a a b b -d b a a b +/3e fde a a d a fe + (d a a b c -d b a a c )e b c + /3e fde a a We bf ) (4.9) 
C = K- 1 p(e def (d a a bd ~d b aa d )(S a e S b f + 28 b f e a e + e a e e b f ) 

+(3(a a e a b f ~aJa b e )(S:S b +2S b f e a e + e a e e b f )) . (4.10) 

The hrst order of these constraints are the linearized constraints of General Relativity. Linearized 
constraints arising from a first class constraint set are Abelian JO]) this applies also to our constraints. 

Now we have to choose our clock variables. Having in mind that we have to invert the matrix 
A K j(a,(j') := {T K (a),Cj(a')} order by order, it would be convenient (although not necessary) if the 
zeroth order of this matrix would be given by the identity matrix; i.e. if the clocks would be exactly 
conjugated to the linearized constraints. 

In metric variables such clock variables were already given in the seminal paper pQ and used in |14j 
to construct the lowest order ground wave function for quantum gravity. Hence we will call these clock 
variables ADM clocks. Their interpretation ^3] is the following: In the gauge where these clock variables 
vanish (or are constant) the metric is in a coordinate system which is as near to the Cartesian one as 
possible. Non-constant clock variables mean that the coordinate system is deformed from the optimal 
one. 

These clock variables can be transformed from the metric variables to the Ashtekar variables. This 
will result in functions which are first order and higher in the fluctuations - we will keep only the first 
order terms. Additionally one has to construct a set of clock variables conjugated to the linearized Gauss 
constraints: 

G T a = (3- 1 e a9e (-W~ 2 d g d d )e d e + w~ 2 (~^W~ 2 d a d d d e + ^d a 8 de )a de (4.11) 

v T a - -W- 2 {~W- 2 d a d d d e + ^d a 5 de ~d e 5 a e ~d d 5 a e )e de (4.12) 
C T = (^W- 2 {-d c t cde e de -p(5 de -W~ 2 d d d e )a de ) (4.13) 
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where W := \/—d a d a . Together with the linearized constraints these clock variables form a set of 
canonically conjugated variables, i.e. { G T a (a), ^G b (a')} = 5%5(a,a') and so on. Here the superscript 
W to the left of a symbol denotes its first order. The Poisson brackets between two clock variables vanish. 

A tensor mode decomposition 4 of the linearized constraints and the clock variables reveals that these 
include all tensor modes except for the symmetric trace-free transverse (STT) modes, which are the 
physical degrees of freedom for linearized gravity: 



(1) G b = 


— 1 / c o (LT r>ac , LL r>ac\ de , n (LT r>dc , TL r>dc , AT r>dc\ a f 

k {6 bc d a ( P de + P de )e +l3e bcd ( P af + P af + P af )a J 


= 




wc = 


K - 1 (2(3e abd d a AT P b f d e a fe ) 


Grpa 


p- l b aS t Sh c LT P h de e, de + W- 2 (S fa d c LL P d f + \5 cb d a T Pt)a de 


V rj-ia 


- W -\-d h LT P b d : - o b tl p£ + \d a T P c / e s cf - \d b LL P b d a e )e de 


C rp 


{Ap)- 1 W- 2 {-d c e eab AT Pte de + (3(-5 cb 7 'p£ - 25 cb LL P%) a de ) . 



Indeed one can introduce a new canonical coordinate system, with coordinates given by the ^TT-modes 
of a ab and e cd , the linearized constraints and the clock variables. We will call these variables STT- modes, 
C-modes and T-modes respectively. 



4.1 Asymptotic conditions 

Since we want to work with asymptotically flat spacetimes we have to formulate asymptotic conditions 
for our phase space variables. These asymptotic conditions have to ensure |151 that the symplectic 
structure is well defined. One such choice is to impose the conditions J^l 

pab 

^ + 0(r- 2 ) (4.15) 

where r is an asymptotic spherical coordinate, B ab is a smooth tensor on spatial infinity with odd parity 
and F ab a smooth tensor on spatial infinity with even parity. 

With this choice of boundary conditions it is well defined to integrate the Gauss constraint with a 
smearing function A J which is 0(r~ 2 ) and the vector and scalar constraints with smearing functions 
which are 0(r _1 ). However we will also have to deal with smearing functions which are 0(r _1 ) and have 
even parity in leading order for the Gauss constraint and are O(r ) with odd parity in leading order for 
the vector and scalar constraint respectively. Additionally we will have the case of a constant smearing 
function for the scalar constraint. Such smearing functions arise if we consider the clock variables as 
smearing functions. In order to make the constraints with such smearing functions functionally differen- 
tiable, i.e. in order to be able to perform integration by parts, we will follow the usual strategy ^3 l*R)| 
and subtract from the constraints those boundary terms that arise if one performs integration by parts. 

Note that there is only a divergence problem for the lowest order part of the constraints, i.e. the 
linearized constraints. This holds also for the various new constraints Ck and Ck ~ the divergence 
problem exists only for the linear order part and these parts coincide with the linear order parts of the 
original constraints. Hence it is sufficient to consider this part of the constraints. Boundary terms which 
arise from the higher order parts of the constraints vanish due to the asymptotic conditions. The smeared 

4 For the definition of tensor modes see appendix 1X1 
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linearized constraints with boundary terms subtracted are 



«G[A] := k- 1 f K\d a e a b + fc bde a ed )da-K- 1 [ A b e\dS a 

is JdT 

= k- 1 j (-d a A b )e\ + A b f3e bde a ed da (4.16) 
Jt 

WV[N] := k- 1 f N a (d a a b b - d h a a b )da - K~ l f (N a a b b - N b a b a )dS a 
Jt JdT 



= -K- 1 ((d a N a )a b b - (d b N a )a a b )da (4.17) 
Jt 

WC[N] := 2n- 1 f3 [ 6 abc Nd a a bc da-2n- 1 f3 [ Ne abc a bc dS a 
Jt JdT 

= -2 K - 1 /3 / e abc {d a N)a bc da (4.18) 
Jt 

where dS a — n a dS, dS is the volume element of the sphere at infinity and n a is the outward pointing 
unit normal and has odd parity. Hence the boundary terms actually vanish for A 6 even and AT, N a odd. 
(Here it is understood that the integration is performed over some coordinate ball with finite radius and 
that one then considers the limit r — > oo). 

The linearized part of C[N] vanishes for a constant N = 1 smearing function. Also we have that for 
such a smearing function the boundary term J__ e abc a bc dS a is equal to the ADM energy modulo terms 
that are proportional to the Gauss constraints 1121 ITf)| and therefore it is a Dirac observable. On the 
constraint hypersurface this ADM energy is equal to ^ 2+ - ) C[l] = J s ^ 2+ ^Cda. 

Henceforth we will understand that we will use the scalar constraints with boundary terms as in 
(|4.18() . Hence we can perform integration by parts if we have a smearing function approaching a constant 
for r — > oo. For smearing functions with stronger fall-off at infinity or for A b even parity and 0(r _1 ), 
N,N a odd parity and O(r ) the boundary terms vanish. 



5 Dynamics 

The clock variables and the expanded constraints from the last section allow us to compute perturbative 
complete observables order by order. First we have to determine the constraints 

C k (<t)~ £do>C j (a>)(A- 1 )t K (a>,a) . (5.1) 

Here the index j runs from to 6, Co := C, C\ to C3 is equal to the vector constraints and C4 to C§ to 
the Gauss constraints. The index K runs also from to 6 with T° := C T, and so on. 
For (|5.1|l we need the inverse of the matrix 

Af{a,a f ):={T K (a),C j (</)} = Sf 6(a, a') + {T K (a), ^C,(a')} 

= : SfS(a,a')+Bf(a,a') (5.2) 

where ^ 2+ ^Cj(a') denotes all terms of order two or higher of the constraint Cj(a'). Hence A^(cr,a') is 
to lowest order given by the (infinite dimensional) identity matrix. This allows us to invert the matrix 
order by order using an iterative equation: 

{A-^ M (a,a') = 5 l M {a,a')- J^ l K Bf (a,a")(A-%(a" ,</)) da" . (5.3) 

Hence we can determine the constraints Ck up to a finite order in the fluctuations in a finite number 
of steps. To calculate the constraints Ck defined in sectional to some order r we have to perform (r — 1) 
times the iteration step described in section [21 For instance the second order of the constraints Ck is 
given by 

^C K {a) = ^C K (a) + [ da' [ da" ™ $ L > (a, a' , a")^ C Ll (a') (1) C i2 (a") (5.4) 
Jt Jt 
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with 

^^ L %a,a\a'')=~{T^(a'),{T^(a''),^C K (a)}} . (5.5) 

This can be used to calculate the complete observable F\f-m with parameter choices t k = to some 
arbitrary finite order m. 

\ F[f , T] {r K =Q) = jr-J da 1 ---da r {---{f^C Kl (a 1 )},---}, [m - r+ ^C Kr (a r )}x 

- — n J S 

(-IfT^^O-.-T^K) . (5.6) 



r=0 



However we are also interested in dynamical questions, that is complete observables for varying clock 
parameters r. Introducing non-vanishing clock parameters into the series for the complete observable 
(13.1(1 we see that it is now a power series in (t k — T K ) which includes the zeroth order term r K . Hence 
the complete observable to the m-th order is not a finite sum anymore. 

Let us consider the power series for complete observables for non-vanishing clock parameters t in 
more detail, separating terms with different powers of T K : 

°° i r 

F [m (r) ~ V-r / d<ri---cM- ••{/,£*>!)},■.• ,CW r (a r )}x 

(r^(^i) - T^iaJ) ■ ■ ■ (r^(a r ) - T^{a r )) 
1 f 

=- H~i / dcr l-" d(7 r{---{/,CV 1 (cri)},--- .(^(ov)} X 

X)^ / d<r 1 -"d<r r {---{/,fi'jf 1 ((ri)},---,(5 J f p (o- r )})x 

r--l ' ^ E 

r 

£ r Kl (ai)T K > (<7 2 ) ■ • ■ (-T^K)) • • • r K "(a r ) + 



9=1 



°° 1 f 

/ d C r 1 ---d < 7 r {---{/,^ 1 (cr 1 )},--- ) CK r ((r r )} 



£ T^aOr^M • • • (-T K <{a q )) ■ ■ ■ {-T Kp (<r p )) ■ ■■r K ^{a r ) 

q=l,p=2,q<p 

+ ... (5.7) 

Since the constraints Cxi (c) Poisson commute up to terms at least quadratic in the constraints among 
themselves, one can rearrange the constraints Ck { (c) in the formula above in an arbitrary order, changing 
the expression only by terms proportional to the constraints. This allows us to write 1(5.7(1 in two ways. 
On the one hand 

F [f , T] (r K ) ~ ]T- / da 1 ---da r {---{f,G Kl (a 1 )},--- ,C Kr (a r )} x 

r=0 T ' ^ S 

r Kl (^i)^ 2 (^)---r^K) + 

X) ? _ 1)| / dfJl ' ' ' d(7 r{- ■ -{fi&Ktivi)},- ■ ■ ,C Kr (<J r )}(x) X 
r=l ^ ^ s 

^(ai)^'^) ^(ov-iX-T^oy)) + 

1 /" 

X 2'fr- 2V / da i--- d < T r{---{f,C Kl (vi)},--- ,C Kr (a r )}{x)x 

r=2 S 



r^(^)^ 2 (^)---T^- 3 (<7 r _ 2 )(-T^-K_i))(-T^( CTr )) 



(5.8) 
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which after a relabelling of the summation index r can be recognized as 

F im(r K )^F [a ^ {fyT] (r K = 0) . (5.9) 

Here tig (/) denotes the evolution of / with respect to the constraints Ck(&) and the parameters T K (a), 
given by the first summand in 1)5. 8[l . 
In the same way one can arrive at 

°° 1 f M 

F[f; T ](T,x) / dax---da q ol t Gm ({•••{/, C Kl (<ri)}, C Kq (v q )}) (x) x 

g=0 q ' ^ S 

{-T K ^a 1 )){-T^{a 2 ))---{-T K ^a q )) . (5.10) 

Consider the complete observable in the form l|5.9|) . There one has first to evolve the phase space 
function / with respect to the constraints Ck(ct) and the parameters t k (a) and afterwards to calculate 
the complete observable. If one is interested in an m-th order approximation, then the difficulty arises 
that in general the m-th order a^, (/) will contain infinitely many terms involving arbitrarily high order 

of the constraints. The reason for this is, that the lowest order in the Poisson bracket y- m 'g, Ck{&)} of 
an m-th order function ( m )<? with the constraints is (to — 1), hence one looses one order due to the first 
order part of the constraints. 

However we are not interested in the complete observable for arbitrary parameter values t k (<t), for 
dynamical questions it is sufficient to be able to calculate complete observables for a one-parameter family 
of parameters. For our case a natural choice is T°(t; a) = t and T K (t\ a) = for all K ^ 0. The parameter 
t € IR would correspond to a notion of time that is as near as possible to the Minkowskian time of the 
background and as we will see later on corresponds to time translation at infinity. For this choice of 
parameters the evolution (/) becomes 

oo , r 

a% o (f) := oT^if) = £ - {■ ■ ■ {/, C [l]}, • • • , C [l]} (5.11) 

r=0 

where C*o[l] := J^C(a)da — 2k~ 1 (3 J as e ahc ab c dS' a according to the definition in section FTTI Now the 
first order part of Cq[1] is equal to the first order part of C[l] and hence vanishes according to equation 

(EH). 

Therefore we are left with an evolution of / with respect to a generating function ( 2+ '(7[l] which is 
at least second order. This allows us to compute the m-th order approximation to this evolution in a 
systematic manner - the highest order term required of ( 2+ '(7[l] is the (to + l)-th order term (for / a 
first order quantity). 

In order to calculate the m-th order approximation of the complete observable (|5.9I) it is sufficient 
to have the m — th order approximation of a^, (/): The complete observable (for parameters t k set to 
zero) of an n-th order quantity is at least of n— th order, that is one does not loose any orders in the 
second step of the calculation in l|5.9() . Moreover the calculation of the m-th order complete observable 
(again for parameters t k set to zero) requires at most (to + 1) terms in the power series for complete 
observables. 

Before considering explicitly the second order approximation of the complete observable l|5.9|) we will 
remark on the formula l|5.10|l 

«£({■■ (*!)},■•• x 

(-T^MX-T^)) • • • (-T Kq (cr q )) (5.12) 

of the complete observable and bring it into a form where the ADM energy appears explicitly as the 
generator for the evolution of the complete observables in the time parameter t, that is as the physical 
Hamiltonian J7j. To this end choose also here the one-parameter family r°(t, er) = t and r K (t,a) = 
for K 7^ 0. The first order term ^C^l] vanishes, that is we can replace oi^' in l|5.12|l by a( 2+)( j ■ Now 

the functions ' 2+ - ) C'o[l] Poisson commute with the clock variables T K (a) since ( 2+ )(7o(<r) does not contain 
any C-modes, i.e. does not depend on the variables canonically conjugated to the clock variables. 



g=0 q ' JE 
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Hence a\ 2+)6 ^T K (a)) = T K (a) and we can write (j5.12j) as 

F um {t)~a\ 2+)0o \F [m {t = ^ , (5.13) 

where we use t as an abbreviation for the one-parameter family r°(t, a) =t and t k (t, a) = for K ^ 0. 
In equation Ij5.13|l we evolve the Dirac observable F^. t T\{t — 0) with respect to the generating function 
( 2+ )Co[l]. However since we evolve a Dirac observable we can add any combination of constraints to 
the generating function, that is we could also use ( 2+ )C[l] which differs from ( 2+ )(7o[l] only by terms 
proportional to the constraints. 

Now ( 2+ )C[l] (and ^ 2+ -'C'o[l]) coincides on the constraint hypersurface with the ADM energy. This 
shows that the ADM energy appears as the physical Hamiltonian, that is as the generating function for 
the time evolution chosen here. Since the ADM energy generates time translations at infinity, we see that 
our choice of time parameter corresponds to time translations at infinity. 



6 The second order approximation 

In this section we will consider explicitly the second order approximation of the complete observable 

F [f . T] {t)~F H (/);T] (Q) (6.1) 

where / is a first order phase space function and we will assume that it commutes with the linearized con- 
straints and with the clock variables. For pure gravity we could choose / = STT P^a a b{u) =: STT a c d(cr) 
or / = STT P^e a b{cr) —■ 5TT e cc i(o'). The higher than second order approximations can be obtained in a 
similar manner. To simplify notation we will introduce H := ( 2+ )C"o[l]. The calculation of (|6.1() proceeds 
in two steps, first one has to calculate the evolution of / with respect to H and afterwards one has to 
compute the complete observable corresponding to a^.(/) with the clock parameters t k set to zero. The 
complete observable -F^t] (0, x) of a phase space function of order to is at least of order to, hence one 
does not loose any orders in the second step. This means that in order to calculate (|6.1[) to second order 
we also need o^(/) to second order: 

co r r 
r=l ' s=0 

= "*«*(/)+ E ^ {{{/.^Ip,^}^}. • (6-2) 

p,q=0 ^ q >' 



Using the identity 



(q+p+l)\ J q\ p\ 
for q,p natural numbers, we can rewrite the sum in (|6.2J) as 



" qw> -rw* (6 . 3) 



""«'*(/) = «!.,*(/)+/ *' £ iJ -4 }lt -{{{J, l2, H}„' 3 'H} m H} r 



p,q=Q 



U 6 (f) + / dt' a<$ s ( {a{l;p(f), ^H} ) . (6.4) 



For higher order calculations one has to use the identity l|6.3|) iteratively. Note that a[ 2)i ^ is the 
propagator for a linear field theory, hence can be given explicitly. 

Finally we have to compute the second order complete observable corresponding to ^Og(f). The 
second order complete observable corresponding to a first order function g has three summands: 

[2] F [g[T] (0) = g+ [ dtnigWClKMK-TKfa)) 

+i f Aa 1 Aa 2 {{g^C Kl {<T 1 )}^C K2 {a 2 )}{-T K -{a 1 )){~T K ^a 2 )) . (6.5) 
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If g is a second order quantity, we only need the first two summands and replace ^Ckx{ cf i) by 
W CKti&i) there. Alternatively since we have that ([3]) 

*W- S2 ;T](0) = F [gi;Ti (0) ■ F [g2]Ti (0) 

%+ 92i T](0) = F [5i:T] (0) + F [g2;T] (0) (6.6) 

we can calculate the second order complete observable of a second order function by computing the first 
order complete observable of its first order factors. 

Note that the complete observables associated to g and g' , i.e. F[ g;T j(0,x) and F[ g ,. T ](0,x) coincide 
weakly if g and g' coincide weakly on the (gauge fixing) hypersurface {T K = 0}. Hence before we 
calculate the complete observable associated to o^(/) we can set all terms vanishing on the gauge fixing 
hypersurface to zero. Therefore we also need to determine a^(/) only modulo such terms. 

6.1 The two lowest orders of the Hamiltonian 

Here we will consider the second and third order of the Hamiltonian H in more detail. We know that 
any finite order of H := ( 2+ )(7o[l] does not contain any C-modes. On the other hand, as explained in 
the last section, we can omit in cn^.(/) any terms containing T-modes. We need therefore to consider H 
only up to terms containing T-modes. Let us denote by H' the function obtained from H by setting the 
T-modes to zero. Then we can use H' instead of H in the propagation of / and furthermore we know 
that the only gravitational modes appearing in any finite order of H' are the STT-modes. 

These consideration simplify very much the calculation of at least the lower order terms of H' since 
one can omit all terms which contain either C-modes or T-modes. This is the advantage provided by 
introducing the constraints Ck as compared to the constraints Ck- 

Consider the second order of H. The second order of Cq[1] is given by 

^C [l] = ^C[l]+ [ Wc j (o')W(A- 1 %(o>, a )do'd* . (6.7) 
For the second order of H = ( 2+ )Co[l] we have to add a term quadratic in the first order constraints: 
(2) C'o[l] = (2) Co[l]-^ |{T^(^i),{T i2 (^), (2) C'o( ( T3)} (1) CL 1 (a 1 )«C L2 (a 2 )da 1 d ( T 2 d ( T3 . (6.8) 

Hence the second order of H is given by the second order of C[l] plus terms which contain C-modes. 
Therefore we can obtain the second order of H' by setting all T-modes and C-modes in the second order 
of C[l] to zero. 

For the third order we need 

( 3 >C [1] = < 3) C[1]+ / WC J (a')W(A- 1 ) (<j' 1 a)da'da + 

/ ^C j {a')^{A- 1 %{a' i a)6a'Aa . (6.9) 

Here the second summand on the right hand side might turn out to be relevant for H' . However as 
shown in appendix 151 the functions J ( - 1 )(A~ 1 )q(<7 / , a)da contain only T- and C-modes. Furthermore 
one can check that the terms one has to add to ^Cofl] to arrive at the third order of Co[l] are at least 
linear in the C-modes. Therefore we can obtain also the third order of H' by restricting ( 3 ^C[1] to the 
STT-modes. However the fourth order of H' will differ from ^ 4 ^C[1] restricted to the STT-modes. In 
particular in this order terms with the inverse derivative operator W~ 2 = (— 9 a 9°) _1 will arise, leading 
to a non-local time generator H' . This reflects the non-locality of our choice of time function. 

7 Gravity coupled to a scalar field 

In this section we will consider gravity coupled to a scalar field and compute the complete observable 
associated to the scalar field to second order. Here the scalar field is assumed to have only small deviations 
from the zero value, that is the scalar field and its conjugated momentum will be counted as phase space 
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functions of first order. As we will see, the first order complete observable coincides with the expression for 
the scalar field on a fixed Minkowski background. Hence we will compute the lowest order gravitational 
correction to this expression. 

In order to couple a scalar matter field to gravity we have to add the following matter contributions 
to the vector and scalar constraints l|4.4|l : 

<V 6 = -wd b <f> (7.1) 

7 

C = i (vr 2 + qq ab d a c^d b <j) + qm 2 ^ 2 ) (7.2) 

where </> and ir denote the scalar field and its conjugated momentum with {(j)(a), Tr(a')} = 7<5(<r, a'). Here 
7 is a coupling constant for the scalar field and m is the mass for the scalar field. We will assume that 
these fields fall off as 0(r~ 2 ) for r — > 00. Furthermore q ab is the inverse metric and q = det(q a b) the 
determinant of the metric. With 

qq ab = !3 2 E*E) q = 3 det(E?) (7.3) 

the scalar matter constraint (|7.2|) can be expanded to 

*C = -!- U 2 + {S ab + 2e ab + e ad e b d ) d a ^d b cf> + 
27 

(1 + e\ + \{e a a e b b - e\e b a ) + ie o6c e de > e ad e be e cf ) m 2 cp 2 ) . (7.4) 

The lowest order of l|7.4[) is the second order and coincides with the Hamiltonian for a scalar field on 
a flat space-time. The third order has a gravitational correction. We need this third order for the third 
order of the propagator H' as defined in section RTT1 As explained there we just need to restrict the third 
order of (|7.4|) to the STT-modes, hence the matter contribution to the third order of H' is 

m<t>H> = lsTT e ab 9a ^ _ (75) 

7 

Let us calculate the first and second order of the propagated scalar field according to (|6.3() . The first 
order is given by 



r=0 



S(t, a; 0, a')Tr(a') + S'(t, a; 0, <j')<f>{cr') da' (7.6) 

where we introduced the propagators 5* and 5" for the scalar field. These propagators and the propagators 
G, G and G" for the gravitons are reviewed in appendix[n] For the second order of the propagated scalar 
field we find 



{2} a%,(4>(a)) = / dt'a ( ( ^,({a[^?(0(a)),( 3 )j?'}) 
Jo 

= / dt' af a) *,( / 2S(t-t',a;Q,a')( STT e ab d a d b ^)(a')da' ) 
Jo Js 

= / dt' [ da' 2 S(t - t', o-; 0, a 1 ) x 
Jo Js 



s 

(G'? d (t', a'; 0, a") STT a cd (a") + G"? d (t', a'; 0, a") STT e cd (a"))da" 



did* \ {S(f, a'; 0, a"')ir(a'") + S'(t\ a'; 0, a'")<P(a'"))da'" . (7.7) 
Js 

The two contributions l|7.6|l and l|7.7|l to the propagation of the scalar field can be summarized by 
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introducing new propagator functions for the scalar field that depend on the gravitational variables 
S gr (t,a;0,a') := S{t, a; 0, a') + 

f At' f da" 2 S(t,a;t' ,a") al' 2)fl ( STT e ab (a"))df S(t' ,a";0,a') 
Jo Js 

S' gr {t,a;0,a') := S(t, a; 0, a') + 

/ At' f da" 2 S(t, a; t\ a") ai m [ STT e a \a"))df df S'(t', a"; 0, a') 
Jo Jr, 

(7.8) 

such that 

^a%{<p{a)) + ^a%{^{a)) = f S gr (t, a; 0, a')ir(a') + S' gr (t,a; 0, a')<j>{a') Aa' . (7.9) 



(2) 



Finally we have to compute the second order complete observable associated to the first and second order 
propagated field. Note however that the second order l|7.7|) is already invariant under the constraints 
modulo terms of second order, since it is a sum of products of two phase space functions which are 
invariant to first order. Hence we only need to compute the second order term corresponding to the first 
order propagated field l|7.6J> . According to ((6.5J) we have 

F[M«W));T\ = ~ I da'(S(t,a;0 7 a')d b (7r v T b )(a') + S'(t,a;0,a')(d b cj ) )(a') v T b (a')) 

da' (Sit, a; 0, a')d a ((d a <f>) c T)(a') - m 2 <j)(a') c T(a') + 

S'(t,a;0,a')7t(a') c T(a')) . (7.10) 

The last term in (|6.5|l vanishes in this case, because in ^'oh, (4>(a)) there only appear matter fields and 
no gravitational fields. 

Now the second order complete observable is given by 

[2] F M „ hT] (t) = V>afe,M*)) + {2) a%Ma)) + {2) F [{1)a%M{a)) , T] (7.11) 

where the explicit expressions for the quantities on the right hand side can be found in equations 
(j7.6l7.7l7.1t)(l . The first term in (|7.11|l coincides with the expression for a scalar field at time i on a 
flat space-time. The other terms contain corrections due to the coupling to gravity: The second term 
^o:*j, (<j>(a)) is due to the fact, that the scalar field propagates on a space time with (non-interacting) 
gravitons, the last term ensures gauge invariance to second order. 

To facilitate the interpretation of the result (|7.11(l note that the first two terms arise also if we evolve 
the scalar field <f> with the time-dependent Hamiltonian 

H{t) = ( da^(ir 2 (a) + (S ab + 2 STT e ab (t))d a q>d b <j>(a) + m 2 2 (a)) (7.12) 

where the time dependence of STT e ab (t) is given by 
s ™e"\t,a) = a\ 2)kl { STT e»\a)) 

(G'? d (t, a; 0, a') STT a cd (a') + G"f d (t, a; 0, a') STT e< d (a')) da' . (7.13) 

Here, if one evolves with the Hamiltonian (|7.12|) one does not treat the gravitational variables as 
dynamical anymore, i.e. the Poisson brackets between the gravitational variables are set to zero. Indeed 
(I7.12|l can be interpreted as the first order (in the graviton field) approximation of the Hamiltonian for 
a scalar field propagating on a graviton background. The first two terms in (|7.11|) are also the zeroth 
and first order approximation to the propagation of the scalar field with the time dependent Hamiltonian 
(|7.12|l . Hence we captured in (|7.11|l the lowest order effect of a scalar field evolving on a graviton 
background. (The last term in (|7.11(l vanishes on the hypersurface, where all gravitational modes except 
for the STT-modes vanish.) Therefore the new propagator functions defined in l|7.8Jl are the zeroth 
and first order approximation to the propagator functions for a scalar field propagating on a graviton 
background. The higher order terms which arise if one evolves the scalar field with the Hamiltonian 
(17.12(1 can be found as a subset of the higher order terms in the perturbative expansion of the complete 
observable associated to the scalar field. 
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8 Poisson brackets: Space— time algebra of observables 

Here we will consider the Poisson brackets between two second order complete observables ^ F^^-t] (ti) 
and m F [<p{a2) . T] (t 2 ). Note that the Po isson bracket of two second order gauge invariants is an invariant 
of first order, we therefore need to consider only zeroth and first order terms of the Poisson bracket. The 
zeroth order will coincide with the result for the field commutator on a flat space-time, in particular 
it will vanish if (tx,ai) and (t 2 ,a 2 ) are space-like related (with respect to the Minkowski metric). The 
first order correction will be a function of the gravitational variables. Fluctuations in the gravitational 
variables will be reflected in fluctuations of the light cones, i.e. the causal structure. 
The zeroth order of the Poisson bracket can be found to be 

= -7 / (S(h, ai; 0, a') S'{t 2l a 2 ; 0, a') - S'^a^O, a 1 ) S(t 2 , a 2 - 0, a'))da' 

JY 

= -jS(ti,ai;t 2 ,a 2 ) ■ (8.1) 

For the first order of the Poisson bracket we have to consider the Poisson brackets between the first 
order term in H7.11fl and the two second order terms in (|7.11|) . To begin with we will show that one of 
the two contributions vanishes: 

{^a%M-i))/ 2) F [Wa%w ^ T] } 

= { (1 >a£,(0(°l)). / {^a%M^)), [2] C K (a')}(-T K (a'))da'} + 



{ / {^a%M^)), l2] CK(<7')}(-T K (a'))da', ^a%(<j>(a 2 ))} 

JY 

/ { {1] a t ?M^)), l2] C K (a')}(-T K (a'))da'} + 

JY 

{ [ {^a%(cb(a 2 ))^C K (a')}(-T K (a'))da', W a %M*i))} 



= (8.2) 

where we used the Jacobi identity and the fact that the matter fields Poisson commute with the clock 
variables in the second equation. 
The other contribution is 

{«4,Wi)), Wa%,(<p(a 2 m + {Wa%Mvi)), W <4M* 2 ))} 

= {Ua%,(<f>(<ri)), f t2 dt' [ 2S(t 2 ,ar,t\a')ag 6 X STT e a \<T'))a^dfai )6 ,(^'))d<T'}- 



o 



((ti,(Ti) «-> (t 2 ,a 2 )) 



= -7/ df / 2S(t 2 ,a 2 ;t',a')a{' 2m ,( STT e ab (a'))d:'dZ' S^a^t' ,a')da' + 

JO JY 

7 I' dt' f 2S(t 1 ,a 1 ;t', C r')a t { ' 2)f{ ,( STT e ab (a'))d: ' d( ' S(t 2 ,a 2 ;t' ,a')da' 
Jo Jy 

= - 7 /' 1 dt / f 2S(t 1 ,a 1 ;t',a')a t ( ' 2)fl ,( STT e ab (a , ))d:'d^'S(t',a , ;t 2 ,a 2 )da' . (8.3) 

J to JY 



Comparing this result with the definition of the propagator function S gr in (|7.8|l we see that for t 2 = 
we can write 

{ l2] F MaihT ](ti), [2l ^(. 2 );T](ii)} =-jS gr (t 1 ,a 1 ;0,a 2 )+O(2) (8.4) 

where 0(2) refers to terms of second order. This result is similar to the Poisson bracket for a scalar 
field on flat space; the flat propagator function is replaced by the "effective" propagator function S gr , 
which to the lowest non-trivial order takes into account the effects of the graviton background. Hence 
we can say that to this order the observables have the causality properties of field observables on such 
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a (fixed) graviton background. Higher order terms of Poisson bracket between two complete observables 
will be much more complicated, for instance already the second order terms will include Poisson brackets 
between the (phase space dependent) propagator functions. Moreover the second order terms may include 
inverse spatial derivative operators making locality considerations quite difficult. 

However we can say, that all the higher order terms vanish in the zero gravity limit, that is for 
k, a a i,, e a b — > 0. Indeed the higher order terms are all proportional to 

7 (grav. field)" 1 (-(matter field) 2 )" (8.5) 

7 

for some powers m, n £ N. We want to compare this behaviour to the case where one chooses four scalar 
fields as clocks for the scalar and vector constraints, see for instance [Sj. To this end we will consider 
the Poisson brackets restricted to the gauge fixing hypersurface, i.e. the hypersurface where the clock 
variables coincide with the parameters r . Furthermore we will consider only infinitesimally separated 
clock parameter values t and t + e. Since we have [1] 

{F [f , T] (r),F [g , T] (r)} = F [{f , g} . : T] (r) (8.6) 

where {•,•}* is the Dirac bracket with respect to the gauge T K (a) = t k and [5] 

F [f . T] (r + e) = F [f+{fJ 6k ^ } ;T] (r) + 0(e 2 ) (8.7) 

we have to examine 

{0(^0, {0(a 2 ), f e K C K da}}* =s M<n), {0K>, / e K 6 K da}} 
Js Je 



+ / {^ 1 ) ! T i (a')}{C' i (a') ! {^ 2 ), / e K C K da}} 



E 

{^Cl^OHTVU^), f e K C K d*}} , (8.8) 



where we used the definition of the Dirac bracket, see for instance [S]. We will be interested in the case 
where the field <p{cr) Poisson commutes with the clock variables, so that the second term in l|8.8|l vanishes. 
Using the property of the Abelianized constraints, that {T L (a), Ck(o'')} — S^6(a, a') + 0(C 2 ), one can 
see that also the third term vanishes: 

{TV), { 0M, / e K C K da}} = {c/>(* 2 ), {T l ((t'), f e K C K da}} ~ . (8.9) 
With the definition (|5.4|l of the constraints Ck 

C K (a) = [ Qia'^A-y^^da' + [ ^ (a, a', a")C Ll (a')C L2 {a")da'da" + 0(C 3 ) 
JS Je 

^ L *(a,a',a") = ~{T^(a'),{T L %a"),C K (a)}} (8.10) 
the first term is equal to 

WK),{^ 2 ),/ e K C K da}} ~ / e K (a){^a 1 ),{^a 2 ),C 3 (<j')}}(A-%(a',a)dada' + 

e K (a){cj ) (a 1 ),(A-%(a\a)} {<f>fa), Cj(a')} dada' + 



E 



e K (a){cl ) (a 2 ),(A- 1 y K (a',a)} i^), C 3 (a')} dada' + 
e K (a){cj ) (a 1 ), C Ll (a')}{4>(a 2 ), C L2 (a")} x 



E 



{T (a'), {T (c"), C i (&'")}} (A~ )k(<j"' , a) dada'da"da'" (8.11) 



18 



We will apply this formula to the case were one uses four scalar fields S (a) as clock variables. We 
have to add the following terms to the scalar and vector constraints: 

1 3 

S C = —y2((Il K ) 2 + qq ab d a S K d b S K + qV(S K )) 
la * — ' 

K=0 

1 3 

s Vb = _ V U K 8 b S K (8.12) 
a * — ' 

where Hk is the conjugated momentum to S K with {S K ,Ul} = aS^ and V(S K ) is a potential for the 
scalar field S . Then the second and third term in (|8.11() vanishe and we are left with 

{^U^), f e K C K da}} ~ 7^1,^2) / e K (a)(A- 1 )° K (a 1 ,a)da + 

a f e if (a){0(a 1 ),C Ll (a')}{0(a 2 ),C L2 (a')} ( 5 LlL2 (A- 1 )O f (a',a)dad ( 7' 

(8.13) 

Here e'(er') := J e K (a)(A^ 1 ) ( j i (a l ,a)da is the translation to the factor that would arise if we would use 
/ e'(a')C(a')da' as time generator. We therefore see 

7^ J eV){^i)A>OH0MA>')}<5 LlL2 d^ ^ e'(a'){^ 1 ),Q(a')}{^a 2 ),C fe ( ( T')}x 

(A-'y Li (a')(A-i)l(a')S L ^da' (8.14) 

as the genuine effect of this choice of clock variables on the Poisson brackets. In l|8.14|l we used that 
the inverse matrix A -1 can be written as (A~ 1 ) J Li (a 1 , a") = (A~ 1 ) 1 Li (a')S(a', a") if one uses scalar 
fields as clock variables. Since A? (a, a') — {T K (a),Cj(<r')} the term (|8.14f) can be interpreted as the 
generalization of the expression "energy of the observed field divided by the energy of the clock variables" 
which we derived for the case of parametrized particles in the introduction, section ^ The correction 
(|8.14(l can be made small by choosing a to be very small (or equivalently the energy of the clock variables 
very large). But one has to keep in mind that through the coupling to gravity backreaction terms arise 
which scale with (positive powers of) n/a, hence one has to balance between the term in (|8.14(l and the 
backreaction terms. In contrast to this result, the corrections to the Poisson bracket in the case of the 
ADM clocks scale in the same way as the backreaction terms. Another point is that using the ADM clocks 
one can perturb around flat space, whereas if one uses scalar fields one has rather to perturb around a 
phase space point with a finite energy density and hence a non-vanishing gravitational field due to the 
constraint equations. 



9 Summary and discussion 

In this work we introduced a perturbation scheme for the calculation of Dirac observables. For this 
perturbation scheme one has to choose a fixed phase space point providing a background space-time. 
Dirac observables can be calculated order by order in the fluctuations around this phase space point. 

We applied this method to general relativity and chose as the fixed phase space point the Minkowski 
background. However the method is also applicable to other backgrounds, for instance a cosmological 
background. Furthermore one can also choose another set of clock variables, as long as these clocks 
provide a good parametrization of the gauge orbits near the fixed phase space point. Our choice of 
the ADM clocks was guided by the aim to obtain observables which approximate very well the field 
observables on a flat background. 

The first order approximation of the complete observables coincides with the observables of the lin- 
earized theory, which gives a precise understanding of how the observables and the dynamics of the 
linearized theory, as for instant the graviton are connected to the observables of the full theory. 

The second order terms of the complete observables associated to matter fields take into account 
the propagation of these matter fields on a graviton background, that is the scattering of matter from 
gravitons, the third order term will contain among other things backreaction terms. 

5 If we would rescale our gravitational variables to a' . := k~ 1 ^ 2 a a f ) and e' ab := K _1 / 2 e a i, the first order of a complete 
observable associated to a matter field is proportional to k° and in general the m-th order to /t( m_1 '/ 2 . 
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Formula l|5.11|l and l|5.13[) provide two different ways to calculate these Dirac observables. In (|5.11|l 
we first time evolve the (gauge variant) field in question and afterwards calculate the gauge invariant 
extension of the result. It is shown in section HTTl that as the generator for the time evolution we can choose 
the function H' , which does only contain the "physical modes" of the linearized theory, i.e. only matter 
fields and the symmetric transverse traceless (STT) modes of the gravitational field. This facilitates the 
interpretation of the resulting time evolution as one involving the dynamics and the scattering of gravitons 
(the STT-modes) . We have different interaction processes contributing to the time evolution and for 
each individual interaction process one could in principle calculate the fully gauge invariant extension. 
In this way one can associate a Dirac observable to each interaction process. 

On the other hand in (|5.13() we first calculate the complete observable or gauge invariant extension 
associated to the field in question and then evolve the resulting observable with a generating function 
which is also a Dirac observable and is given by the ADM energy. This shows that our time evolution 
corresponds to time translations at infinity and that with our choice of time we have a Hamiltonian 
whose matter part coincides with the matter part of the Hamiltonian constraint. Note however that 
this time generating function is not unique - one can add an arbitrary combination of the constraints 
without changing the result on the constraint hypersurface. For instance one can choose between the 
ADM energy, which in our case is of finite order and the function H which includes infinitely high order 
terms. 

The resulting Dirac observables have a local interpretation stemming from the flat space-time limit, 
they give the fields at some space time point coordinatized by (t, a). The interpretation of the coordinates 
(t, a) is not as straightforward as in the case where one uses matter fields as clock variables: what one 
can say is that the measurement of these observables has to be performed with respect to a reference 
system which is as near to the Minkowskian one as possible. Thus to zeroth order the proper distance 
between two points can be obtained by using the flat metric, for higher orders one has to use the (time 
evolved) gravitational field e ab (t, a). However a better understanding of the geometrical meaning of these 
coordinates would be helpful. An interesting question is whether the ADM clocks can be understood as 
an approximation to a set of clocks which has a more obvious geometrical interpretation. 

Furthermore we suggested to consider the Poisson algebra of the complete observables corresponding 
to different space-time points, in order to learn more about the locality properties and the interpretation 
of the complete observables. We made some preliminary steps into this direction and calculated the first 
order of the Poisson bracket. The first order approximation of the Poisson bracket l|8.4(l can be interpreted 
to reflect the causality properties of a space-time with (non-interacting) gravitons. Hence to this order 
one can say that the complete observables with this choice of clocks represent a local measurement. 
But higher orders will include inverse derivative operators stemming from our choice of clock variables. 
In particular our choice of time parameter t is very global since it corresponds to time translations at 
asymptotic infinity. It may help to define another one-parameter family of clock value parameters t k (f), 
for which the t k change only in a finite region, effectively introducing a boundary in space-time. 

Choosing scalar fields as clock variables will lead to a better locality behaviour, since the complete 
observables associated to these clocks will not involve inverse derivatives, see also jHj for a discussion of 
the advantages using scalar fields. However we showed that the Poisson brackets between (infinitesimally 
time-separated) complete observables and hence the uncertainty relations between these complete ob- 
servables include a term which is inverse to the energy of the clock fields. If one wants to keep this term 
small one has to increase the energy of the clock fields which leads to larger backreactions and hence to 
a bound from below for the additional term in the Poisson bracket. 

In summary we think that the approximation scheme introduced here can be very useful to calculate 
Dirac observables and to explore their properties. A better understanding of the properties of complete 
observables is needed, in particular of the question how the choice of clock variables influences these 
properties. 

A Tensor modes 

Similar to the longitudinal and transversal modes for a vector field on K 3 one can introduce tensor modes 
for a tensor field. For a proof of the completeness of these modes, see ^B]- To begin with we define the 
projector onto the transversal modes of a vector field by 

(p.v) a :=p b a -v b :=(S b a + W- 2 ■ d a d b ) ■ v b . (A.l) 
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This allows us to introduce the following basis of tensor modes: 

( LT P ■ T) ab = (S c a - p c a ) ■ v d h ■ T cd 2 left long, right transv. modes 

{ LL P ■ T) ab = (S c a - pi) ■ (St - pf) ■ T cd 1 left and right long, mode 

[ TL P ■ T) ab = p c a ■ {Si - pi) ■ T cd 2 left transv. right long, modes 

( T P ■ T) a b — \pab ■ p cd ■ T c d 1 symm. transv. trace part mode 

( AT P ■ T) ab = \ (p c a ■ pi - P c b ■ pt) ■ T cd 1 antisymm. transv. mode 

{STTp , T ^ b= i . p d + p c . p d _ pabp cdj . 2 symm. transv. tracefree modes (A.2) 

Using the projector property p ■ p = p, it is easy to see that the projectors x P are orthogonal to each 
other and satisfy P ■ Y P = S XY x P. Furthermore the set of projectors is complete, that is 

J2 Xp z£ = sc Ji ■ ( A -3) 

x 

B First order of the matrix A 

Here we display the relevant terms to find the third order of the time generator H' used in section 16.11 

{ G T a (a), WC[1]} = W~ 2 /3{-2e aef d e d 9 + 3e 9ef d e d a )a fg 
{ v T a {a), WC[1]} = 2/3((W- 4 d a d c d e + W- 2 d c Si + W- 2 d e Si)a ce + 

/3W- 2 (d d d b e abc - d a d b e dbc )e dc ) 

{ c T{a),^C[l}} = - l -W~ 2 d a d e e a e + l -W- 2 d c d c e\-pW- 2 e efd dea fd (B.l) 

The right hand sides of the equations l|lj.l|l do not contain STT-modes, which shows that the third order 
of H' cn be obtained by just restricting the third order of C[V[ to the 5TT-modes. 

C Propagators 

Here we will review the propagators for the scalar field and the graviton on a flat background. We will 
start with the scalar field. Given initial values 4>(u) and ir{a) on E the scalar field <f>(t, a) at a later time 
t can be calculated to be 

cj>(t,a) = jr^^Cll}}/- 

■r=0 

= cos [(-A CT + m 2 fl 2 t\ 0(a) + (-A CT + m 2 )- 1 ' 2 sin [(-A CT + m 2 ) 1 ' 2 ^) (C.l) 

where ^ 2 ^C[1] = L i(7r 2 + 9 Q 09 a + m 2 (j) 2 ) and A = S ah d a d b . Introducing a delta function this can be 
rewritten as 

4>{t,(j) = / cos[(-A CT +m 2 ) 1 / 2 i] ( 5((7,CT')</'(^) + (-^+™ 2 )" 1/2 sin[(-A CT + m 2 ) 1 / 2 ^(a, ( T')^(^)dff' 

JT, 

S'{t,(j;0,a')(j){<T') + S(t,<T;0,a')ir(<T)d<j' . (C.2) 

s 

From the last equation one can read off the propagator functions S and S'. This definition of the 
propagator function can be generalized by 

S(ti,ai;t 2 ,a 2 ) := 5(*i - t 2 , a\] 0, a 2 ) (C.3) 

and the analogous definition for S' . We have that S(ti, a±;t 2 , a 2 ) = — S(t 2 , a 2 ; t\, <j\) is odd under 
permutation of t%,t 2 , whereas S'(ti,ai;t 2 ,a 2 ) = S'{t 2 ,a 2 ;ti,ai) is even. 

In the same way one can determine the propagator functions for the graviton field. The Poisson 
brackets of the graviton fields STT e a b = STT Pa b e cd and STT a ab with the second order Hamiltonian 

^H' = - [ {2(5e bed STT e\ d b STT a ad - (i 2 STT a ab STT a ab ) da (C.4) 
« Jt. 
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defined in section 16.11 are given by 

{ STT a ab (a), ^H>} = 2/3D% STT a fe (a) := 2p\d c e cd %5 db 8{ + 6 d j() STT a fe {a) 
{ STT e ab (a), ^H'} = 2/3 2 STT a ab (a) - 2(1D% STT e fe {a) 
{ STT e ab (a),^H'} 2 = [{2PD)XI STT e/e(a) . (C.5) 

Hence the graviton fields are evolved according to 

" TT a{t,a) = f G f a t(t,a;0,<j') STT a fe (a')da' 
Is 



STT e(t, a) = J G'il{t, o; 0, a') STT a fe (a>) + G" f a t(t, a; 0, a') STT e fe (a')da' (C.6) 



with 

G f a e b (t,a;0,a') = [expp/JiA,)]^,^) 

fe. 



00 f2r+l 

G'it(t,a;0,a') = 2f3 2 ]T — — — [(2/3A,) a *-]&*(a, a') 
G" f a e b (t,a;0,a') = [ cxp(-2/ft£> ff )] %8{a, a') . (C.7) 
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